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POLYNOMIAL-TIME  IDENTIFICATION  OF  OPTIMAL  ROBUST  NETWORK 
FLOWS  UNDER  UNCERTAIN  ARC  FAILURES 


VLADIMIR  BOGINSKI  AND  CLAYTON  W.  C'OMMANDhR 


Abstract  Network  flow  problems  have  a  wide  variety  of  important  applications  in 
many  areas.  Although  deterministic  formulations  of  these  problems  are  well-studied,  in 
many  practical  situations  one  has  to  deal  with  uncertainties  associated  with  possible  fail¬ 
ures  of  network  components  (e.g.,  each  arc  has  a  probability  of  failure).  Formulations  and 
optimal  solutions  of  these  problems  need  to  be  adjusted  to  take  into  account  these  uncer¬ 
tainty  factors.  The  main  difficulty  arising  in  addressing  these  issues  is  the  dramatic  increase 
in  the  computational  complexity  of  the  resulting  optimization  problems.  We  propose  LP- 
based  solution  methods  for  network  flow  problems  under  a  set  of  failure  scenarios,  which 
allows  for  robust  solutions  to  be  found  in  polynomial  time.  We  justify  this  fact  by  proving 
that  for  network  flow  problems  under  uncertainty  with  linear  loss  functions,  the  number  of 
scenarios  required  to  approximate  the  mean  of  the  loss  distribution  for  any  fixed  e  >  0  with 
probability  (1  -  a),  for  a  C  (0. 1],  is  polynomial  with  respect  to  the  size  of  the  network. 


1.  Introduction 

Network  flow  problems  are  among  the  most  well-studied  topics  in  operations  research 
with  numerous  papers  and  textbooks  devoted  to  their  study  [1,  3,  4,  9,  14,  15,  17,  20], 
These  problems  have  a  wide  variety  of  important  applications  in  many  areas.  Although 
deterministic  formulations  of  these  problems  are  well-studied,  in  many  practical  situations 
one  has  to  deal  with  uncertainties  associated  with  possible  failures  of  network  components 
(e.g.,  some  arcs  have  a  probability  of  failure).  Formulations  and  optimal  solutions  of  these 
problems  need  to  be  adjusted  to  take  into  account  these  uncertainty  factors. 

The  main  difficulty  arising  in  addressing  these  issues  is  the  dramatic  increase  in  the 
computational  complexity  of  the  resulting  optimization  problems.  For  example,  Corea  and 
Kulkami  [6]  have  considered  the  minimum  cost  flow  problem  where  the  length  (cost) 
of  each  edge  is  a  random  variable.  They  construct  continuous  time  Markov  chains  to  derive 
stable  algorithms  computing  the  distribution  of  the  minimum  cost.  The  stochasticity  of  the 
edge  weights  makes  the  problem  particularly  challenging  in  that  the  size  of  the  state  space 
grows  exponentially  with  the  size  of  the  graph.  Glockner  and  Nemhauser  [12J  have  also 
considered  dynamic  network  flow  problems  where  the  arc  capacities  are  random  variables. 
They  describe  novel  decomposition  techniques  for  the  resulting  multi-stage  stochastic  lin¬ 
ear  program.  Doulliez  and  Rao  [  1 0]  have  studied  the  problem  of  determining  the  maximum 
flow  in  a  graph  subject  to  a  single  arc  failure.  Along  this  line,  Aneja  et  al.  [2]  considered 
the  maximum  flow  PROBLEM  and  present  a  strongly  polynomial  algorithm  for  maxi¬ 
mizing  the  residual  flow  in  a  graph  after  a  single  arc  is  destroyed.  Cormican  et  al.  [7] 
have  formulated  stochastic  programming  formulations  of  the  network  interdiction 
PROBLEM  in  which  arcs  are  removed  in  order  to  minimize  the  expected  maximum  flow  on 
the  network.  These  problems  are  modeled  as  two-stage  stochastic  programming  problems 
and  are  fAf.f'-hard. 
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Figure  1 .  A  graphical  depiction  of  VaR  and  CVaR  [5]. 

In  this  paper,  we  propose  linear  programming  based  solution  methods  for  network  flow 
problems  under  uncertainty.  In  particular,  we  provide  a  polynomially  solvable  formulation 
for  the  ROBUST  minimum  COST  Fi.ow  problem  (rmCF),  in  which  each  arc  has  a  prob¬ 
ability  of  failure.  To  our  knowledge,  this  is  the  first  formulation  of  its  type.  Moreover, 
we  prove  for  linear  programming  problems  under  uncertainty  with  linear  loss  functions, 
the  required  number  of  scenarios  needed  in  order  to  ensure  that  the  sample  mean  is  within 
e  >  0  of  the  true  mean  with  probability  (1  a),  for  0  <  a  <  1 ,  is  polynomial  in  the  size  of 

the  input.  To  our  knowledge,  this  is  the  first  time  such  questions  have  been  considered  for 
optimization  problems  under  uncertainty. 

2.  Scenario-Based  Approaches  for  Optimization  Problems  under 

Uncertainty 

Optimization  problems  that  involve  uncertainties  are  typically  formulated  using  random 
samples  (or,  scenarios)  from  loss  distributions  [8].  An  important  question  that  needs  to  be 
answered  is,  how  many  scenarios  are  needed  to  accurately  reflect  the  statistical  properties 
of  the  loss  distributions  with  a  high  confidence  level ?  More  specifically,  does  the  required 
number  of  scenarios  grow  polynomially  or  exponentially  with  respect  to  the  size  of  the 
input?  Clearly  this  is  an  important  problem,  and  to  the  best  of  our  knowledge  has  not  been 
previously  addressed  in  the  literature  on  optimization  under  uncertainty.  In  this  paper,  we 
will  formulate  the  robust  minimum  cost  flow  problem  under  uncertainty  and  prove 
that  for  network  optimization  problems  with  linear  loss  functions,  for  any  fixed  positive 
e  >  0  and  a  g  (0,  Ij,  O  (nr )  scenarios  arc  sufficient  to  ensure  that  the  sample  mean  of 
the  total  loss  associated  with  possible  edge  failures  is  within  e  from  the  true  mean  with 
probability  (1  —  a)  (where  m  is  the  number  of  edges  in  the  network). 

Next,  we  discuss  the  basic  statistical  concepts  (sometimes  referred  to  as  risk  measures ), 
which  we  utilize  in  this  study. 

2. 1 .  Quantitative  Loss  Measures.  One  of  the  most  well-known  risk  measures  used  in 
robust  optimization  under  uncertainty  is  known  as  Value-at-Risk  (VaR)  [  1 3  ] .  VaR  provides 
an  upper  bound,  or  percentile  on  a  given  loss  distribution.  For  example,  consider  an  appli¬ 
cation  in  which  a  constraint  must  be  satisfied  within  a  specific  confidence  level  a  6  (0, 1], 
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Then  the  corresponding  a- VaR  value  is  the  lowest  value  £  such  that  with  probability  a, 
the  loss  does  not  exceed  ^  [16].  In  economic  terms,  VaR  is  simply  the  maximum  amount 
at  risk  to  be  lost  from  an  investment.  VaR  is  the  most  widely  applied  risk  measure  in 
probabilistic  settings  primarily  because  it  is  conceptually  simple  and  easy  to  incorporate 
into  a  mathematical  model  [5].  However  with  this  ease  of  use  come  several  complicating 
factors.  Some  disadvantages  are  that  the  inclusion  of  VaR  constraints  increases  the  num¬ 
ber  of  discrete  variables  in  a  problem.  Thus  a  polynomially  solvable  problem  is  likely  to 
become  hard  [11]  after  the  VaR  constraints  are  added  to  the  model.  Also,  VaR  is  not 
a  so-called  coherent  risk  measure,  implying  among  other  things  that  it  is  non-convex  and 
not  sub-additive. 

Another  risk  measure  closely  related  to  VaR  is  the  so-called  Conditional  Value-at-Risk 
(CVaR).  CVaR  is  a  more  conservative  measure  of  risk,  defined  as  the  conditional  expec¬ 
tation  of  the  loss  under  the  condition  that  VaR  is  exceeded.  Rockafellar  and  Uryasev  [18] 
proved  several  important  results  regarding  optimization  of  CVaR,  which  make  this  risk 
measure  rather  attractive  from  the  optimization  viewpoint.  In  particular,  CVaR  has  been 
shown  to  possess  the  properties  that  VaR  lacks;  in  particular,  it  is  coherent  (which  includes 
convexity  among  other  properties).  This  makes  this  statistical  measure  much  more  con¬ 
venient  to  handle  in  optimization  models.  A  graphical  representation  of  the  relationship 
between  CVaR  and  VaR  is  shown  in  Figure  1.  In  order  to  define  CVaR  and  Var  we  need 
to  determine  the  cumulative  distribution  function  for  a  given  decision  vector  subject  to 
some  uncertainties.  Suppose  L(x,y)  is  a  loss  function  associated  with  a  decision  vector 
x  €  X  C  K",  and  a  random  vector  y  £  R"'  which  is  the  uncertainties  that  may  affect  the 
performance.  Assume  that  y  is  governed  by  a  probability  measure  P  on  a  Borel  set,  say  Y 
[5].  Then  the  loss  L(x,y)  for  each  xeXisa  random  variable  having  a  distribution  in  R 
induced  by  that  of  y.  Therefore  the  probability  of  L(x,y )  not  exceeding  some  value  C,  is 
defined  as 

y(*,Q  :-P{y\L(x,y)<Q.  (1) 

By  fixing  x,  the  cumulative  distribution  function  of  the  loss  associated  with  the  decision  x 
is  thus  given  by  y(x,  Q  [21]. 

Given  the  loss  random  variable  L(x,y)  and  any  a  e  (0,1),  we  can  use  equation  (1)  to 
define  a- VaR  as 

CaM  :=min{^€R:y(x,Q  >a}.  (2) 

From  this  we  see  that  the  probability  that  the  loss  L(x,y)  exceeds  £a(x)  is  1  -  a.  Using 
the  definition  above,  CVaR  is  the  conditional  expectation  that  the  loss  according  to  the 
decision  vector  x  dominates  Ca(x)  [18].  Thus  we  have  a-C’VaR  denoted  as  <j>a(x)  defined 
as 

<Mx)  :=  E{L{x,y)\L(x,y)  >  £«(*)}•  (3) 

In  order  to  include  CVaR  and  VaR  constraints  in  optimization  models  we  must  charac¬ 
terize  £a(x)  and  <|)a(x)  in  terms  of  a  function  Fa  :  X  x  R.  •— »  R  defined  by 

Fa(x,Q  :=C+  (max  (Z.(x,y)  -  £,0}}.  (4) 

In  [18],  Rockafellar  and  Uryasev  prove  that  as  a  function  of  Fa(x,  Q  is  convex  and  con¬ 
tinuously  differentiable.  Moreover,  they  show  that  a-CVaR  of  the  loss  associated  with 
any  x  6  X,  i.e.,  <j)a(x),  is  equal  to  the  global  minimum  of  Fa(x,Q,  over  all  £  6  R.  Fur¬ 
ther,  if  Aa(x) argminrtR  Fa(x,  Q  is  the  set  consisting  of  the  values  of  C,  for  which  F  is 
minimized,  then  Aa(x)  is  a  non-empty,  closed  and  bounded  interval  and  i[„(x)  is  the  left 
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endpoint  of /lu(*)-  In  particular,  it  is  always  the  case  that  £a(x)  S  argminr(  K  Fa(x,Q  and 

VaW  ~  Fa(x,ia(x))  [18]. 

This  result  gives  a  linear  optimization  algorithm  for  computing  a-CVaR.  It  is  a  result  of 
the  convexity  of  Fa(x.Q,  that  we  are  able  to  minimize  CVaR  for  x  6  X  without  having  to 
numerically  calculate  §a(x)  for  every  x.  This  has  been  show  by  Rockafellar  and  Uryasev 
in  [19].  Further,  it  has  been  shown  in  [19]  that  for  any  probability  threshold  a  and  loss 
tolerance  C,  that  constraining  4>o-.  (*)  <  C  is  equivalent  to  constraining  Fa(x,  Q  <  C. 

3.  Minimum  Cost  Flow  Problem  under  Uncertainty  with  CVaR 

Constraints 

An  instance  of  the  minimum  COST  flow  problem  (MCP')  consists  of  a  directed  graph 
G  =  {V.E),  where  each  edge  (arc)  (ij)  €  E  has  an  associated  cost  cy  per  unit  of  flow  along 
this  edge,  as  well  as  a  capacity  u  -,j  denoting  the  maximum  amount  of  flow  that  can  traverse 
edge  (/,  /).  For  each  node  i  €  V,  d,  denotes  demand  (supply)  of  node  i.  Then  the  MCE  can 
be  formulated  as  the  following  linear  program 

(MCF)  min  £  cy*y  (5) 

('V)efc' 
s.t. 

X  XU  ~  X  xj‘  = v ' e  K-  (6) 

0<j Cij  <Ujj,V  i,j  FV.  (7) 

Suppose  that  L(x,y)  is  a  loss  function,  where  w{0, 1 }  is  a  vector  of  decision 

variables.  In  the  context  of  a  network  flow  problem,  x /,■  represents  the  total  amount  of  flow 
through  arc  Further  let  y  be  a  random  vector  representing  the  edge  failures  in  the 
graph  defined  as  follows: 


Vij  ■- 


1 ,  with  probability  ptj, 

0,  with  probability  1  -  p,j, 


where  p,}  is  the  probability  of  failure  for  the  arc  (/,_/). 
Furthermore,  for  each  scenario  s  =  1 , 2, . . . ,  S, 


yh  '■= 


1 ,  with  probability  pih 
0,  with  probability  1  —  p,j. 


Then,  the  random  variable  C(x,_v)  representing  the  total  loss  is  given  as 

n 

L{x,y)  -  Y  Wo- 
ij- 1 

Further,  for  each  scenario  s  —  1 , ...  ,S, 


E'(x.y)  =  £  Xij/ij. 

‘j- 1 

Then  the  characteristic  function  Fa(x,  Q  for  CVaR  constraints  is  given  as 

Ea(xf)  =  C  +  7 ~  X  7C'  max{  X  X  xvrfj  ~  C.o), 

1  Us_-1  li=lij=l  J 


(8) 


(9) 


(10) 


(ID 
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where  7tv  is  the  probability  of  scenario  y1,  for  s  —  1 , 2, . . . ,  S.  Since  L(x,y)  is  linear  with 
respect  to  x ,  F„( x,  Q  is  convex  and  piecewise  linear  f  1 6].  If  we  assume  that  each  scenario 
is  equally  likely,  that  is  tts  =  5,  V s  =  1,2, ...  ,S,  then  (11)  reduces  to 


Fa{x'Q  ~ c+ t  max{  .t  -  &0}-  (12) 

Then  we  can  model  the  the  robust  MCF  problem  with  CVaR  constraints  as 


(RMCF-1)  Minimize  cyxy  +  ^  M,pt  (13) 

{(<V)ef)  ieF 

s.t. 

X  X>J~  X  xfi  =  d<  -  A'V  i  6  v,  (14) 

{J:{iJ)eE}  0':l/V)€£} 

0  <  Xij  <  Ujj  V  (ij)  6  £. 

C  +  (1  -a)5^maX{  (!5) 

£gR,  (16) 

Pi  >  0,  V  /  €  K.  (17) 


Note  that  the  mass-balance  constraints  in  (6)  from  formulation  MCF  have  been  relaxed 
in  formulation  RMCF-1.  This  modification  is  made  in  order  to  ensure  that  arc  failures 
do  not  render  the  instance  infeasible.  Notice  that  penalty  terms  representing  the  possible 
loss  of  flow  (e.g.,  unsatisfied  demand)  at  node  i  have  been  added  to  the  objective  function 
and  the  corresponding  constraints.  Clearly,  it  is  desired  that  the  penalty  variables  p,  are 
equal  to  zero  in  the  optimal  solution,  which  is  modeled  by  multiplying  these  variables  by 
sufficiently  large  coefficients  M,  in  the  objective. 

We  can  linearize  Fa(x, Q  by  using  dummy  variables  ts,  s  —  1,2,... ,5,  and  replacing 
Fa(x.  Q  by  the  linear  function  £  +  Ej=i  h  and  adding  the  set  of  linear  constraints 


is  >  Vj  =  1,2,. ,.,S, 

ts  >  0,  V.9-  1,2,. ..,5. 


(18) 

(19) 
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Thus,  the  fully  linearized  mathematical  programming  formulation  of  the  robust  min¬ 
imum  COST  FLOW  problem  with  Conditional  Value-at-Risk  constraints  is  given  as 


(RMCF-LP)  Minimize  Cijxy  4-  £  A/,/?,  (20) 

( i,j)cF.  iev 
s.t. 

X  x‘j  -  X  xji  ~  di  Pi,VieV,  (21) 

<22) 

C,  V  J=  1,2,. ..,5,  (23) 

U=  I 

ts>0,  Vs=l,2,...,S,  (24) 

c  e  R,  (25) 

p,  >  0,  V  i  6  V,  (26) 

0  <  JCy  <  u/y,  V  f,y  6  E.  (27) 


Note  that  this  LP  formulation  contains  0(\V\  +  |£j  +  S)  variables  and  0(\V\  +  |£|  +S) 
constraints.  Therefore,  the  subject  of  further  investigation  is  how  large  the  number  of 
scenarios  S  needs  to  be  to  ensure  sufficiently  good  statistical  properties  of  the  optimal 
robust  solutions.  In  the  next  section,  we  address  this  issue  and  identify  the  conditions 
under  which  S  is  polynomial  with  respect  to  \E\  (and  |F|),  which  makes  the  proposed 
problem  formulation  polynomially  solvable. 


4.  On  thf.  Required  Number  of  Scenarios 


In  this  section,  we  discuss  the  issues  of  obtaining  a  “good"  approximation  of  the  true 
distribution  of  the  losses  in  the  above  formulations.  More  specifically,  we  investigate  the 
issue  of  the  required  number  of  scenarios  that  would  ensure  that  the  true  mean  of  the  loss 
distribution  is  close  enough  to  the  sample  (scenario-based)  mean  with  a  high  confidence 
level.  In  particular,  the  subject  of  special  interest  would  be  the  dependency  of  the  required 
number  of  scenarios  on  the  input  size,  e.g.,  the  size  of  the  network.  This  is  particularly 
important  as  the  number  of  scenarios  considered  will  directly  affect  the  number  of  addi¬ 
tional  constraints  in  the  model.  In  particular,  if  the  required  number  of  scenarios  needed 
is  exponential  in  the  size  of  the  input,  then  the  polynomial  solvability  of  the  problem  is 
undermined  by  the  exponential  number  of  constraints.  In  this  section,  we  prove  that  the 
required  number  of  samples  needed  to  accurately  approximate  the  mean  of  the  distribution 
of  a  linear  loss  function  is  in  fact  polynomial  in  the  size  of  the  input  (e.g.,  the  number  of 
nodes  and  edges  in  the  network)  under  certain  conditions. 

Let  the  sample  mean  loss  be  given  as 


L\  +  Li  H - YLs  1  ,  *  v'  V-  ,t 

L  ■= - c - =  =  cX  X 

S—  1  i-1  !J— 1 


(28) 


Also,  let  sample  mean  of  each  component  of  the  random  vector  y,j  be  defined  as  yij 
5  Zf-i  y]j-  Let  the  true  mean  of  the  loss  distribution  be  ul-  With  this,  we  have  the  following 
proposition. 


Proposition  1.  Consider  a  linear  loss  function  for  a  network  flow  problem,  such  as  the  one 
defined  above.  Then  for  any  e,  a  G  R  such  that  £  >  0,  and  a  6  (0,1],  the  required  number 
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of scenarios  S  to  guarantee  that  P{\L  —  ///,|  <  e)  >  1  —  a  is  0(m2/  £2),  where  m  —  |£|  is  the 
number  of  edges  in  the  original  network.  Furthermore,  for  any fixed  £,  S  —  0  (m2). 


Proof. 


P(\L-g,.\<Z)  - 


p(ll2  2  xuy]j  -  2  xupu\  < 

\  J  s-l  i,j~l  !.j-\  / 


p(l  «7')  I  <  e 


s— 1 


(29) 

(30) 

(31) 


=  P[\'L  xu(yu~Fj)\  <ej- 

Now,  let  U  max,y{«y}  be  the  maximum  capacity  of  any  single  arc  in  the  network.  Then 
continuing  from  (3 1 )  we  have  that 

p(\  X  -/»</)!  <eJ  >  P^mUma.x{\yu- inf)  <ej  (32) 

=  />^max{^;-^y|}<^j  >  1-a.  (33) 

Let  i*  and  /’  be  the  indices  for  which  | —  ^,y|  is  maximum,  and  let  oj.j,  be  the  variance 
of  v,-y* .  Then,  using  the  notation  c«  representing  the  critical  value  of  the  normal  distribu¬ 
tion  corresponding  to  the  confidence  level  ( 1  —  a)  (e.g.,  for  1  —  a  =  0.95,  ca  —  1 .96)  the 
minimum  required  number  of  scenarios  S,  needed  to  satisfy  the  condition  above,  can  be 
found  as  follows. 


CaGJ.  j. 

Vs 


e 

mU 

ciaj. U1  ,  ,  , 

- 2 - m2  =  0(m2/V). 


Furthermore,  for  any  fixed  e,  S  —  0  {m2).  Thus  we  have  the  desired  result. 


(34) 

(35) 
□ 


Note  that  without  the  loss  of  generality,  the  result  of  Proposition  1  holds  for  any  loss 
functions  that  are  linear  with  respect  to  the  uncertain  parameters  y,  where  these  parameters 
can  denote  any  uncertainties  in  the  model,  such  as  costs,  capacities,  or  other  application 
specific  parameters. 

The  proof  of  Proposition  1  provides  the  justification  for  the  use  of  Conditional  Value-at- 
Risk  constraints  for  network  flow  problems.  That  is,  by  incorporating  CVaR  constraints  we 
are  able  to  guarantee  that  the  sample  mean  of  the  loss  distribution  can  be  made  arbitrarily 
accurate  while  the  size  of  the  problem  remains  polynomial  with  respect  to  the  input.  Thus 
we  can  find  accurate,  robust  solutions  to  these  problems  in  polynomial  time  using  standard 
linear  programming  techniques. 


5.  Numerical  Example 

In  this  section,  we  present  a  numerical  example  demonstrating  the  effectiveness  of  the 
proposed  approach.  Consider  the  three  graphs  shown  in  Figure  2.  Suppose  we  are  consid¬ 
ering  the  minimum  COST  flow  problkm.  Subfigure  2(a)  represents  the  original  network 
and  the  triplets  (cjj,Ujj,pjj),  correspond  to  the  cost,  maximum  capacity,  and  probability  of 
failure  for  each  arc  (/,_/). 
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(b)  The  ideal  optimal  solution.  The  objective  value 
is  295. 


(c)  The  optimal  solution  with  CVaR.  In  this  exam¬ 
ple,  5  -  100  and  the  objective  value  is  340.2. 


FlGURK  2.  The  min-cost  flow  problem  is  solved  for  the  graph  in  (a). 
Cij,Ujj,  and  p-,j  represent  the  cost,  maximum  capacity,  and  probability 
of  failure  for  arc  (i,j)  respectively. 


In  Subfigure  2(b),  the  optimal  ideal  solution  is  shown.  The  ideal  case  represents  the 
deterministic  MCF  problem,  or  equivalently,  the  probabilistic  version  when  all  edge  failure 
probabilities  are  0.  The  optimal  solution  for  this  instance  is  295  units  of  flow.  However 
by  taking  into  account  the  edge  failure  probabilities,  we  see  that  this  is  not  the  preferred 
solution.  In  2(c),  the  optimal  solution  is  shown  with  CVaR  constraints.  In  this  example, 
the  total  loss  in  the  worst  10%  of  scenarios  is  constrained  to  be  less  than  0.23  units  of 
flow.  We  generated  100  scenarios  uniformly  at  random  for  this  example,  with  up  to  20% 
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of  the  arcs  failing  in  any  scenario.  These  are  arbitrary  values  and  can  be  specified  by  the 
user  as  the  situation  calls.  In  this  case,  the  objective  function  value  is  340.2.  We  see  that 
in  fact,  the  solution  tends  to  push  more  flow  over  arcs  with  lower  probabilities  of  failure 
accepting  the  higher  cost  for  the  hedging.  This  is  exactly  the  behavior  we  would  predict. 
An  interpretation  is  that  by  diversifying  the  flow  across  various  arcs,  we  are  able  to  ensure 
that  the  solution  is  robust  in  the  event  of  (multiple)  arc  failures.  Moreover,  note  that  the 
optimal  solutions  for  robust  formulations  of  the  minimum  COST  FLOW  problem  can  be 
fractional.  Recall  that  for  the  deterministic  version  of  the  mcf  problem,  the  optimal  solu¬ 
tion  values  are  integral  if  all  the  parameters  of  the  original  model  are  integer.  Depending  on 
specific  applications,  one  can  either  operate  with  fractional  flows,  or  use  certain  rounding 
techniques  to  convert  the  flows  to  integer  values. 

6.  Conclusions 

In  this  paper  we  proposed  polynomial-time  linear  programming  based  solution  ap¬ 
proaches  for  network  flow  problems  subject  to  uncertain  arc  failures.  Further,  we  have 
shown  that  for  linear  loss  functions,  a  polynomial  number  of  scenarios  is  sufficient  in  or¬ 
der  to  approximate  the  true  mean  of  the  loss  distribution  arbitrarily  close.  A  numerical 
example  was  given  demonstrating  the  effectiveness  of  the  proposed  procedure  in  the  con¬ 
text  of  the  minimum  cost  flow  problem. 
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